Identifying the Bose glass phase 
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Introducing disorder into the Bose-Hubbard model at integer fillings leads to a Bose glass phase, 
along with the Mott insulator and superfluid phases 1}. We suggest a new order parameter: det(p), 
where p is the one body density matrix, which is nonzero only within the Mott-insulator phase. 
Alongside the superfluid fraction f 3 , it is then possible to distinguish the three phases. The Bose 
glass phase is the only phase which has vanishing det(p) and superfluid fraction. The vanishing of 
det p in the Bose glass phase occurs due to the partial fragmentation of the condensate into localized 
fragments, each with zero superfluid response, which implies the presence of unoccupied sites and 
hence the presence of lines of zeros in the one body density matrix. In the superfluid phase, det p 
vanish for another reason - due to the macroscopic occupation of a single particle state. Finally, we 
suggest the enhancement of the three body decay rate in the Bose glass phase, as an experimental 
indicator for the presence of localized fragments. 
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I. INTRODUCTION 

The onset of superfluidity in random media and its 
connection to the existence of a condensate fraction has 
been an intriguing research subject in the 4 He community 
for many years 013 I n recent years, ultra-cold atomic 
states in optical lattices have provided an experimental 
means to answering these long-standing problems £| . 

The remarkable tunability and dynamical control of 
optical potentials offers the opportunity to emulate 
elaborate physical situations and create novel quantum 
many-body states of fundamental interest. In particular, 
the study of disordered systems comprises an important 
avenue of research. Specifically, the Anderson localiza- 
tion transition, perhaps the best known example for a 
quantum phase transition from an extended to a local- 
ized state, has never been unambiguously observed in an 
experiment. Interacting bosons in a disordered poten- 
tial can present interesting and sometimes unexpected 
behavior, as we shall discuss shortly. It is known 0, [j| 
that interacting lattice bosons in the presence of a dis- 
ordered on-site potential exhibit a phase known as Bose 
glass, characterized as an insulating phase with no gap 
in the excitation spectrum. 

The Bose-Hubbard Hamiltonian p| with on-site disor- 
der is given by 



H 



J 



-h.c- 



U 



i i i 

(1) 

where a,i is the annihilation operator for an atom at site 
i, J is the tunnelling strength, U is the repulsive on site 
interaction strength, and is the random on site po- 
tential, drawn from a flat distribution centered around 
zero, of width A. 

When Wi = for all i, this Hamiltonian exhibits a 
second order quantum phase transition from an energy 
gapped Mott-insulator to a gapless superfluid state at in- 
teger filling (number of atoms equal the number of sites), 
as the ratio U/J is changed The critical value for this 
transition depends on the dimensionality of the underly- 



ing lattice. 

This model was extensively studied by a variety of 
methods, including quantum monte-carlo in ID and 2D 
0,0, exact solutions in ID for small lattices @, renor- 
malization group [jj, and recently in the context of cold 
atoms in optical lattices using the mean field approxima- 
tion [lj} ■ One of the main conclusions of these studies 
was that due to the disorder, there is a new incompress- 
ible yet insulating phase, named the Bose glass phase. 
It was conjectured by Fisher et. al. and verified by 
[Tl| and ^2 that any transition from the Mott-insulator 
phase to the superfluid phase in the presence of disorder, 
should follow through the Bose glass phase (at integer 
filling). 

In this paper we suggest a new criterion for identifying 
the Bose glass phase, based on the comparison between 
the superfluid fraction, and the structure of the one body 
density function. We show that in the presence of inter- 
actions, however small, a localized condensate becomes 
unstable, and fragmentation occurs. In the thermody- 
namic limit and a finite interaction strength, the number 
of localized eigenfunctions scales with the size of the sys- 
tem, so the occupation of each fragment (cigenfunction 
of the one body density matrix) tends to zero. As we 
will explain below, one can use two order parameters, 
the superfluid fraction / s , and det p, where p is the one 
body density matrix, to distinguish between the three 
phases, since f s is non-zero only in the superfluid phase 
and det(p) is non-zero only in the Mott insulator phase 
(where f s = 0). 



II. INTERPLAY BETWEEN KINETIC, 
POTENTIAL AND INTERACTION ENERGY 

In this section we will use a variational argument, to 
demonstrate how the relative strength of the kinetic po- 
tential and interaction energies, can stabilize a different 
ground state, either a Mott insulator (strong interaction 
limit), Superfluid (both disorder and interaction compa- 
rable to the kinetic energy), and Bose glass phase (dis- 
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FIG. 1: Superfluidity as a function of the disorder strength 
for different lattice sizes at zero interaction energy. Note that 
this superfluid fraction is the response of the many-body wave 
function to a phase twist in the boundary conditions, hence 
it is equal to 100% for zero disorder and interaction strength, 
without contradicting the landau definition, which is a stabil- 
ity criterion 

order or interaction or both are dominant). Surprisingly, 
we will find regimes where by either increasing the inter- 
action strength U or the disorder strength A we can in- 
duce superfluidity in the system. The latter phenomenon 
deserves a special name and is called disorder induced 
order (see e.g. IT^V In what follows, we devise a 
variational argument that compare three typical limiting 
forms of the many body ground state, and show how as 
a function of U and A the nature of the ground state 
changes in accord with the known phase diagram of the 
system. We further point out how this result supports 
our suggested criteria. 

Consider first, a II?, non interacting, single particle 
disordered system with the following Hamiltonian: 

H = -_(^ata i+1 + a % a\ +1 ) + W t a\a l . (2) 

i 

As before, Wi is drawn randomly from a flat distribution 
of width A centered around zero. All the eigenstates fa 
of this Hamiltonian arc localized, and we denote the low- 
est energy as Eq and the corresponding eigenstate to be 
</>o, so Hfa ~ Eo4>q. Inserting into this system N non 
interacting bosons, all of them condense in the lowest en- 
ergy state. Trivially, the total many body wave function 
is then given by: ipAG&i, x N ) = Hi <t>o{%i), where AG 
stands for Anderson glass. This state is an example of a 
completely localized state, with zero superfluid response, 
but nevertheless, it is macroscopically occupied. Some 
authors refer to this state, as an Anderson glass [Tfillr^ |. 

In Fig. 1, we plot the superfluid fraction of this 
ground state, as a function of the disorder strength for 
L = 100, 200 and 400 sites. We calculated the superfluid 
fraction f s exactly using the following formula |l§| 

/--w<-5*i*w-£^>. < 3 » 



where T is the kinetic energy operator and J is the cur- 
rent operator. In the following section we elaborate fur- 
ther on this definition of superfluidity. 

It is evident from Fig. 1 that the critical point A c 
where superfluidity disapear tends to A c = in the 
thermodynamic limit L — > oo. This is of course ex- 
pected, since any disorder is strongly localizing in ID 
|l9j | . What happens if we now turn on the interaction to 
a small, positive value? Assuming the usual S-wave inter- 
action, the interaction energy of the localized state fa is 
UoN(N— 1)/2, where Uq = gj dx \fa\ A , and q is the effec- 
tive ID interaction coefficient (2^. This energy is much 
larger than the energy of a fragmented state of the form 
i>BG = Ili fai x i) nz where the integers rii are zero above 
some k < N, and ^ rii — N (as we will see shortly, the 
case k = N corresponds to the Mott insulator phase), and 
fa(x) are the eigenstates of the free particle Hamiltonian 
in Eq. [21 The interaction energy of this state is ap- 
proximately given by J2i U^n^m - l)/2 + J2ijtj UijUiUj, 
where Ui — g J dx \ fa | is the inter state interaction en- 
ergy, and Uij = q J dx\fa\ 2 \faj\ 2 is the intra state inter- 
action energy|2j|. It is not difficult to find a set of n^s 
which gives a smaller energy than a macroscopically oc- 
cupied condensate wave function of the form tJjag- To 
do so, we start by putting the first particle in the lowest 
energy state fa we then add the next particle, searching 
among all states for the state that will minimize the to- 
tal energy (including interaction) . We repeat this process 
until we fill the lattice with N — L particles, to obtain a 
many-body state that has the following form: 

iI)bg(xi,...,x n ) = 1 [[(l) i {x i ) fii . (4) 

i 

Hence, by this method we obtain a specific set of occupa- 
tion numbers n^'s, some of which are zero. We calculated 
the cj)iS for a lattice with 500 sites, which is effectively 
thermodynamic in the regime of A's and J7's that we 
employ. 

Next, we construct two more many body wave func- 
tions: i/)sf(si, xn) = L~ N / 2 which corresponds to 
the superfluid state in the absence of disorder and inter- 
actions, since the condensate wave function then is 1/ \/L, 
and all N particles occupy this extended state. The sec- 
ond wave function is if)Mi{x\, xn) = IL e ( x ~ x i)i 
where e(x — Xi) is 1 at site i and zero otherwise. This 
state approaches the Mott insulator (MI) ground state 
for large U's and no disorder. 

In Fig. 2 we plot the total energy of each of these four 
many boy wave functions: ipAG, ipBG, 4>SF and tpMi as a 
function of the interaction strength U, for A/J = 2. The 
total energy of each of the following states is given by: 
Ebg = J2i n i E i + Uim{ni - l)/2 + UijUiUj, where 
Ei is the eigenstate of the free single particle Hamiltonian 
Hq, Esf = U S FN(N-l)/2-2(N-l)J+Tr{H ), is the 
single body extended state, L is the length of the lattice, 
Usf = 9/L, and H is the free particle Hamiltonian, 
E M i = Tr(H ) - 2N J, 
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FIG. 2: Energies per particle (in units of J), of the four vari- 
ational many body states: ipsF, i^bg, tpMi and t/>ag, as a 
function of the interaction strength U/J, for A/J = 2. The 
ij)AG wave function has an energy which is larger then all the 
other states, for all non-zero A's. From the graph it follows 
that the minimal state is ordered ipBG,ipSF,ipBG,i>Mi, as ex- 
pected. This shows that interactions tend to delocalize at 
small A's, since they penalize localization to specific sites on 
the one hand, and the gain in kinetic energy of being de- 
localized is high, on the other hand. For larger interaction 
strength, this is no longer true, and a superfluid is formed. 
At even larger U's the Bose glass phase is re-entered. Fi- 
nally, when the interactions are even stronger, the particles 
will minimize the energy just localizing at the lattice sites, 
disregarding the disordered potential. 



and Bag = NE +gf dx |0 O | 4 N(N - l)/2. 

We find that there are 4 regimes: for A/J < 1.3 
The minimal energy state corresponds to ipBG- For 
1.3 < A/J < 2.35 we find that the minimal energy 
state is tpSF, which corresponds to the superfluid phase. 
As expected, in this regime, the condensate fraction is 
large, and the condensate wave function is extended. For 
2.35 < A/J < 7.1 the tpBG phase is again minimal. 
Finally, for A/J > 7.1 the state ipMi has minimal en- 
ergy. Comparing these approximate results to the known 
phase diagram we conclude that this simple variational 
argument already exhibit the rich interplay between the 
kinetic energy, the disordered potential energy and the 
interaction energy by displaying the right ordering of the 
phases. Starting from the Bose glass phase and increasing 
the interactions, the superfluid phase can become favor- 
able. For larger interaction strength, the Bose glass phase 
once again stabilizes, and for even larger interactions, fi- 
nally, the Mott insulator phase wins. This correctly de- 
picts the true behavior of the ground state [l^. Also, 
we note that the localized condensate is unstable for any 
non-zero interaction strength, hence, the Anderson glass 
in practice, is a phase which strictly speaking can only 
exist on the U = line. Indeed, when interactions are 



present, the energy cost of concentrating all the particles 
inside the localization length, grows without bounds in 
the thermodynamic limit. 

We note that the one body density function p derived 
from ipSF has one macroscopic eigenfunction, but since 
Tr(p) = N, it follows that some of the eigenstates are 
zero (since the size of p is N by N), this is also true 
when the macroscopic occupation of the extended state 
is not total due to quantum depletion. The structure of 
p for the ipBG state is different, since instead of having 
long range order, the density matrix is composed of small 
blocks, separated by lines of zeros, indicating that there 
are lattice sites which are not occupied, due to their large 
potential energy. In both cases, the determinant of the 
density matrix is zero. Deep in the Mott insulator phase, 
the determinant is one, since the one body density func- 
tion is the identity matrix. Near the transition, still in the 
Mott insulator phase, the ground state wave function is 
approximately |1, 1, 1) + (J/U) J2i ^d\ +1 di\l, 1, 1) 
which also has det(p) > 0, e, < 1 and are equal only 
when A = 0. We therefore assert that det(p) is non zero 
in all the MI phase. This observation is just another 
way to express the fact that only in the Mott insulator 
phase, a gap is formed. When the number of occupied 
eigenstates (fragments) of p equals the number of parti- 
cles, this gap will be proportional to U, which is just the 
cost of moving a particle to an occupied fragment. 



III. DEFINITION OF CONDENSATE AND 
SUPERFLUID FRACTIONS 

Following [23| , we define the condensate fraction as 
the properly normalized largest eigenvalue of the one 
body density matrix p(x, x'). For example, when a single 
particle wave function x( x ) exists such that p(x,x') — 
Nx(x')*x(x), then the condensate fraction is N - the 
number of particles in the system. Note that this defini- 
tion does not imply any specific characterization of the 
condensate wave function x( x )i which can be extended 
or localized. 

Superfluidity is manifest whenever there is a difference 
in the momentum density response to a transverse vs. a 
longitudinal imposed velocity field I23L I24I l2o| . Here, we 
employ the phase twist definition [181 126| , where we im- 
pose a phase twist 9 along one of the dimensions of the 
system, corresponding to the application of a transverse 
velocity field [24]. We equate the superfluid velocity v s 
with — t, where m is the effective mass of the atom, 
and L is the length of the system 24]. If the phase 
twist is small, then the difference in the ground state 
energy between the twisted (Eg) and untwisted (Eq) sys- 
tems can be attributed to the kinetic energy of the su- 
perfluid. The superfluid fraction f a is therefore given 
by (l/2)mv 2 J s N = Eg - E , with N the total number of 
atoms. For a normal fluid, this difference is zero, indicat- 
ing an identical response to longitudinal and transverse 
imposed velocity fields. 
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IV. STUDYING THE EXACT SOLUTION 

Using an exact solution for 8 particles in 8 sites, we 
show that with this criterion we can roughly obtain (up 
to our finite size resolution) the correct phase diagram of 
the system. 



A. A lattice with one defect 

The effect of a random potential on an interacting 
bosonic system can be identified by measuring the su- 
perfluid fraction. In practice, for a finite size system 
with a longitudinal extension L Zl residual superfluidity 
is expected to persist as long as £; oc > L z . Remember- 
ing that in ID any disorder is localizing, we first consider 
the simplest kind of disorder, namely, a single site defect, 
by adding a negative external potential at one particular 
site, with strength A. In Fig. 3, we show the result of 
an exact diagonalization of Eq. Q we performed, for a 
ID system with 8 atoms and 8 sites, and a varying defect 
strength A (in units of J) . From the figure it is apparent 
that starting from the superfluid regime (U/J= 1) and 
increasing A, the superfluid response decreases rapidly 
to zero, while the condensate fraction increases to 100%. 
As we lower the defect, atoms localize there, without the 
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FIG. 3: Exact calculation of the largest eigenstate of p (filled 
circles) and the superfluid fraction (hollow circles) for a one 
dimensional lattice with 8 atoms and 8 lattice sites, and a 
single site defect potential —A. (a) In the weakly interacting 
case (U/J — 1), increasing A/ J causes all the atoms to enter 
the defect, so the maximal eigenvalue of the one body density 
is 100% of the particles, while the superfluid fraction decrease 
rapidly to zero, as the condensate wave function becomes lo- 
calized, (b) Interactions are stronger (U/J = 15), so at first, 
lowering the defect causes derealization of the atoms, hence, 
superfluidity becomes possible. Only at higher values of A/ J 
it is energetically favorable to become localized at the defect, 
and superfluidity disappears. 

ability to participate in the superfluid flow. When all 
the atoms are captured by the defect, they all occupy 
the same single particle wave function which is localized 
at the defect site. Starting, on the other hand, from 
the Mott insulator phase (U/J — 15), and no disorder 
(A = 0), essentially all the atoms are pinned to the lat- 
tice sites. Increasing A causes an onset of the super- 



fluid fraction, since in the presence of large interactions 
atoms can lower their energy by being in an extended 
state, hence making superfluidity possible. This will be 
the basis for understanding the more elaborate statement 
that one can onset superfluidity by increasing the disor- 
der strength - disorder induced order. Only for larger 
A, the potential energy can compensate for the gain in 
interaction and kinetic energies by localizing again at the 
defect site, hence losing superfluidity as indicated in the 
figure for A/ J > 10. We confirm that the condensate 
wave function is indeed extended (comparable non-zero 
values throughout the lattice) whenever superfluidity is 
present pHE^. 

In this section we refereed to the condensate fraction, 
since for a single impurity which is properly scaled it can 
remain finite even at the thermodynamic limit, whereas 
(as will be shown in the next section) for a disordered 
potential the condensate fraction is non zero only for a 
small (mesoscopic) system, and tends to zero in the ther- 
modynamic limit. 



B. Exact solution with a random on-site potential 

Next, we return to the disordered potential, which acts 
on all 8 sites. We exactly diagonalize Eq. (JIJ and cal- 
culate the determinant of the one body density matrix 
and the superfluid fraction, for different values of U /J 
and A/J, where now, again, A is the width of the flat 
distribution from which we draw the random values of 
the on-site potential. The results are shown in Fig. 5. 
Looking at the superfluid fraction as a function of A for 
U = 0, we see that there is a rapid decrease of superfluid- 
ity towards zero due to the strong localization (the super- 
fluidity starts at a non-zero value due to the finite size of 
our lattice). Looking at the dependence on U for A = 0, 
on the other hand, we observe the usual superfluid to 
Mott-insulator transition, which is slightly smeared due 
to finite size. It is also evident from Fig. 5 that det p is 
non zero only at the Mott insulator phase. 

There are two somewhat surprising features apparent 
in Fig. 5. First, starting from small A/J, and increasing 
the interactions, we observe the onset of superfluidity, 
indicating that in this regime, interactions tend to delo- 
calize. Second, starting from the Mott-insulator phase 
(e.g at U/J — 8), and increasing the disorder strength, 
induces superfluidity, a phenomenon called disorder in- 
duced order [l3T |. 

By inspecting the eigenvalues of p we verified that in 
the regimes where the superfluid fraction f s is zero, the 
condensate is fragmented, i.e. the number of non-zero 
eigenvalues of p is larger than one. While both the dclo- 
calizing role of U and disorder induced order are already 
present in the single defect model discussed above (see 
e.g. Fig. 3), fragmentation can not occur for single site 
"disorder" since there is always an overlap between any 
extended state and a localized state. 

Comparing these results to the known phase diagram 



|12| . we find that, whenever the Bose glass phase is 
present, at least one but less than / eigenvalues of p 
(I being the total number of sites), are non-zero while 
the superfluid fraction is zero. In the Mott insulator, we 
observe complete fragmentation, and all the eigenvalues 
are non-zero (with f s = 0). Finally, in the superfluid 
phase, some of the eigenvalues are zero, while the large 
eigenvalue is large (of the order of 1). 

Further insight is gained by looking at specific eigen- 
functions of p. To avoid the limitations of finite size, we 
inspect the Bose glass and the Mott insulator phases in 
a regime in the A/ J — U / J plane, where the localization 
length is < 2a which is less than the length of the sys- 
tem / = 8a, where a is the lattice spacing. In Fig. 4 we 
plot the eigenfunctions of p for the Mott insulator phase 
(MI), The superfluid phase, and the Bose glass phase. 
We note that In the superfluid phase (U/J = A/J = 1), 
the largest eigenvalue of p corresponds to 95% of the 
particles, and the accompanying eigenstate is extended 
(comparable non zero values anywhere on the lattice). 
The last eigenvalue of p are zero i.e. they are unoccu- 
pied. In the Mott insulator state (U/J = 15, A/J = 1), 
we see that all the eigenvalues are non-zero (the largest 
eigenvalue equals 22% of the total number of particles, 
the second largest is 17%, and the smallest is 7%), in- 
dicating complete fragmentation. Also note that the 
eigenfunctions are all extended. Finally, In the Bose 
glass regime, all the eigenfunctions are localized and 
some of the eigenvalues of p are zero, indicating par- 
tial fragmentation. In particular, due to the finite size, 
when the disorder is strong and the interactions are weak 
(A/ J = 15,U/J= 1), we find that the largest eigenvalue 
is 97% and the eigenfunction is localized, with £ ~ 2a. 
When both the interactions and the disorder are strong 
(U/J — 15, A/J = 60), partial fragmentation occurs. 
The largest eigenvalues of p is then 58%, the second, 
third and fourth largest are respectively: 17%, 13% and 
12% and the rest are zero. The results presented here are 
of a singe realization of the disorder, which we choose as 
representative from a larger set of disorder realizations, 
all with similar quantitative results. 



V. A NEW ORDER PARAMETER FOR THE 
MOTT INSULATOR PHASE 

A single, base invariant, order parameter that is non- 
zero only in the Mott insulator phase is det(p). Since 
Tr(p) — N, it follows that only when there is exactly 
one particle in each eigenfunction, all the eigenvalues of 
p are equal to 1 then det(p) = 1, which happens in the 
limit U — > oo. In the superfluid phase and in the Bose 
glass phase, det(p) = 0. We therefore propose that the 
combination of the superfluid fraction and det(p) can dis- 
tinguish the three different phases. 

Inspecting det(p) as a function of U/J for A = 0, we 
find that the superfluid to Mott insulator transition oc- 
curs at 4.6 in accordance with We also verified for 
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FIG. 4: Absolute value squared of the eigenfunctions rji(x) 
of the one body density matrix /3y = (g.s.\a\a,j\g.s.) , as a 
function of position (in units of lattice spacing). Ai is the 
corresponding eigenstate of r/i, indicated only for the largest 
(i — 1) and smallest (i == 8) eigenvalues. 



large U det(p) tends to unity (for example at U — 1000, 
det(p) > 0.9999). 

Thus, the Bose glass phase can be thought of as a phase 
composed of localized eigenstates of p, without phase 
relations between them. The total many body wave 
function is therefore, to leading order, ip(x\, x n ) ~ 
Yli 4>i(xi) ni , where ^rii — N and some of the m's are 
larger than 1. 

In the presence of interactions, a single localized con- 
densate at the lowest point in the potential is unstable 
and will tend to be fragmented to microscopic fragments 
(for large A/J) or to delocalize (for small A/J). This is 
because the interaction energy is proportional to the den- 
sity squared, and thus whenever interactions are present, 
the ground state wave function will take a form that has 
the lowest possible density. A single extended condensate 




FIG. 5: Superfluid fraction, and det(p) as a function of U/J 
and A/J. The Bose glass phase is characterized by zero super- 
fluidity and det(p) = 0, The Mott insulator phase is the phase 
with non zero det p and zero superfluidity, and the superfluid 
phase naturally has non zero superfluidity and det p — 
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is formed (hence invoking superfluidity), when U /J is less 
than the homogeneous Mott insulator-superfluid transi- 
tion (so the energy penalty for occupying more than one 
atom per site, is small compared to the kinetic energy 
gain) and also A/ J is small (so the penalty for occupy- 
ing higher peaks in the disordered potential is also small) . 
Alternatively, fragmenting to a set of non-overlapping lo- 
calized condensates, each centered at a different local- 
ization center, with negligible energy differences (in the 
thermodynamic limit), is favorable when A/J is large. 

Hence we expect that in the Bose glass phase, the 
number of fragmented condensates iVj oc will generally 
be larger than one. Moreover, Ni oc should scale like L d 
where d is the dimension of the underlying lattice. In the 
Bose glass phase Ni oc < N, since if Ni oc = N, i.e. the 
number of fragments equals the number of sites, then 
an energy gap will form, whereas the Bose glass phase 
is not gapped. When the number of localization cen- 
ters is a finite fraction of the total number of sites - a 
mesoscopic system, the number of atoms per localized 
condensate will also be a finite fraction from the total, 
so the condensate fraction will be non-zero. Moreover, in 
a finite system, we expect a regime of parameters (arbi- 
trary disorder and weak interactions) , in which the finite 
gap in the excitation spectrum (that exists due to the 
finite size) will protect a single localized condensate from 
fragmenting. This reg ime was called by several authors 
an Anderson glass |l6LfTi| . We stress that in the thermo- 
dynamic limit, the Anderson glass only exist along the 
line U = (zero interactions). 

Due to the limited number of lattice sites and atom 
number in our exact calculation, in Fig. 5, and contrary 
to the variational calculation presented in section two, 
we could not resolve the Bose glass phase intervening 
between the Mott insulator and the superfluid phases 
P, 0, 0| . Keeping this limitation in mind, we were 
able to recover the correct structure of the phase dia- 
gram, as compared to previous work done with a different 
tool [I3 ■ Our results can also be scaled to larger num- 
bers. To demonstrate this, we repeated the calculations 
for N — 6 atoms and sites, and obtained a qualitatively 
similar phase diagram to that with N — 8. Moreover, 
using finite size scaling in ID in the limit of A = 0, we 
recover the Mott insulator transition at U c — 4.85 J, in 
good agreement with the value 4. 65 J presented in 
We also checked the half filling case for 6 atoms and 12 
lattice sites, where since N < I, det(p) = indicating the 
absence of the MI phase as expected. Finally, we note 
that the fact that in the presence of disorder a Bose glass 
will always intervene between the Mott insulator and su- 
perfluid phases, seems more evident using our criterion, 
since due to the breakdown of translational invariance, 
the only way from the superfluid phase, where a macro- 
scopic condensate exists (no fragmentation), to the Mott 
insulator state, where all the eigenfunction of p are oc- 
cupied (complete fragmentation), goes through gradual 
partial fragmentation. 
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FIG. 6: Exact calculation of the total three body decay rate 
as a function of A/J for different values of the U/J (pre- 
sented in units of (U/J) c ), for N = I = 8 particles and sites. 
Starting at the superfluid phase and increasing A, localiza- 
tion increase the density and hence the three body decay rate 
sharply increases. When U/J is large, the decay rate is still 
notably larger in the Bose glass phase compared to the Mott- 
insulator and superfluid phases, but it is less than the decay 
rate for smaller values of U, due to larger fragmentation that 
decreases the density. When all the particles are localized at 
the lowest defect, the total decay rate reaches its maximal, 
asymptotic value, which here is small due to small size of the 
system (compare with Fig. 7. 



VI. USING THE THREE BODY LOSS TO 
IDENTIFY THE SUPERFLUID TO BOSE GLASS 
TRANSITION 

In this section, we propose a new experimental indica- 
tor for the presence of localized states within the Bose 
glass phase, when interactions are small. To do so, we 
propose that for a mesoscopic system such as the sys- 
tems described i n l28t l2flj . it is possible to use the three 
body decay rate [30, |3l| > to identify the presence of lo- 
calized states. We consider a collection of ID trapped 
BEC's created by loading a 3-D BEC into a deep two- 
dimensional optical lattice, and introduce a controlled 
disorder of strength A relative to Jjioi a more acces- 
sible experimental scheme see [25^, H^)- The ID three 
body loss rate coefficient is r = T 3D /127r 2 er 4 , where 
T 3D = 5.810 _42 m 6 /sec is the corresponding 3D three 
body loss rate coefficient for 87 Kb 31], ay = 85nm 
is the assumed radial confinement, and L = 150pm 
is the condensate length [2^. Hence, we obtain that 
T 1D = 9AlQ- 16 m 2 /sec. In Fig. 6, we plot the total 
three body decay rate T to t as a function of the disor- 
der strength A relative to J, for different values of U/J, 
according to 

Ttot = ^2 {g-s]a\a\a\'aia t a t \g.s.) / L 3 . (5) 

i 

Looking at Fig. 7 together with Fig. 5, we observe 
the following features. When the system is in the Mott 
insulator phase, the total three body decay rate T to t is 
negligible, since the amplitude for three atoms to occupy 



7 




A/J 



FIG. 7: Mean field calculation of the total three body de- 
cay rate as a function of A/J for different values of the UjJ 
(presented in the corresponding (U/J) c units), for 101 atoms 
and sites. Similar to the exact result of Fig. 6, we observe 
that when the interactions are weak, the increase in the three 
body loss rate as a function of the disorder strength A/J is 
pronounced. The decay rate decreases, when the interactions 
increase, since the presence of interaction cause the density 
to decrease due to fragmentation. 

the same site is effectively zero. In the superfluid phase, 
the bosons are now delocalized across the lattice. As a 
result the amplitude to occupy a single site with three or 
more atoms increases to a small value, and T to t also in- 
creases to a small value. The most dramatic effect is seen 
when we enter the Bose glass phase and the interactions 
are weak. In this regime, due to the the existence of lo- 
calized condensates with increased density, we observe a 
notable increase in T to t ■ When the number of fragmented 
condensates increases, we see a corresponding decrease in 
the loss, as expected. The step-like behavior of T tot is due 
to the fact that the capture of another particle is a no- 
table event, since the total number of particles is small. 

In Fig. 7 we repeat the analysis, using a novel mean 
field simulation based on the Gutzwillcr approximation 
[3^, 13 , that we developed, which is capable of working 
at a constant filling factor for arbitrary external poten- 
tials. With this simulation, we also calculated an upper 
bound for the superfluid density, and recovered qualita- 
tively the upper part of Fig. 5. In Fig. 7, we calculated 
the total decay rate from a ID lattice with 101 atoms 
and sites. We calculated the density profile in the pres- 
ence of disorder, and the approximate ground state of the 
system, which we insert to Eq. [5] to obtain the loss rate. 
We note that although the simple Gutzwiller ansatz is 
unable, by construction, to recreate fragmentation, the 
density, profile is expected to be accurate for our pur- 
poses. The results, shown in Fig. 7, are consistent with 



the exact results, up to a trivial scaling of U /J and A/J 
due to the change in the critical point in the mean field 
model. They support our claim that three body loss is a 
good measure for localization in a mesoscopic system, in 
the regime where interactions are weak. 

Finally, we note that localized condensates may also 
be identified by careful examination of their momentum 
distribution observed in time-of-flight images ,8]. This 
can be considered as a complementary measure, to the 
three body loss, since as interactions increase, fragmen- 
tation becomes larger hence the loss rate decreases, 
while the interference is washed out. 



VII. CONCLUSION AND OUTLOOK 

In this paper we propose a new criterion for identify- 
ing the Bose glass phase, as a phase composed of local- 
ized fragmented eigenfunctions of p (with some of the lo- 
calized eigenfunctions remaining unoccupied), and zero 
superfluidity. Starting from the Bose glass phase and 
increasing the interaction strength can cither induce a 
transition to the superfluid phase, where there is one 
macroscopically occupied extended condensate, or to in- 
duce further fragmentation. When the number of frag- 
ments equals the number of sites, we enter the Mott in- 
sulator phase, which also has zero superfluidity, but with 
a non-zero energy gap. A summary of these observation 
is encapsulated by the second order parameter det p that 
we propose, alongside the superfluid fraction f s , which is 
non zero only in the Mott insulator phase. We showed 
that det(p) is non zero in the Mott insulator side, both 
near and far from the transition to the Bose glass phase. 
We also argued that we can capture some of the features 
of the transition, in particular, that the Bose glass phase 
is reentered as U is increased, before entering the Mott 
insulator phase, by using a variational ansatz, which also 
obeys our criterion. We find for an exact solution for a 
small system with 8 particles and sites, that using our 
suggested criterion we can obtain the correct phase dia- 
gram (up to our limited finite size resolution), and that 
fragmentation mechanism occurs for large C/'s. 

Finally, we discuss the possibility to observe localiza- 
tion in a realistic experimental setup in the regime of 
small interaction strength, using the increase in the three 
body decay rate. Several intriguing effects resulting from 
the interplay between disorder and interactions such as 
disorder induced order, arise already in a single defect 
scheme, which may be more accessible experimentally. 

This work was supported in part by the Israel Science 
Foundation, and by the Minerva Foundation. 
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